Anisotropic, Fermi-Surface-Induced Variation in T c in MgE>2 Alloys 
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Using coherent-potential to describe disorder, Gaspari-Gyorffy approach to evaluate electron- 
phonon coupling, and Allen-Dynes equation to calculate T c , we show that in Mg\- X M X E>2 (M = 
Al, Li or Zn) alloys (i) the way T c changes depends on the location of the added/modified k-resolved 
states on the Fermi surface and (ii) the variation of T c as a function of concentration is dictated by 
the B p DOS. In addition, using full-potential calculations for MgMB^, we show that (i) at x = 0.5 
a superstructure can form in Mg\- x Al x B% but not in Mg\- X H X B2 or Mg\- x Zn x B2, and (ii) B 
layer shifts towards the impurity layer, more for Al than for Li or Zn. 

PACS numbers: 74.25.Jb, 74.70.Ad 
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Since the discovery of superconductivity in MgB 2 
the experimental Jl|-^| and theoretical efforts 
have greatly improved our understanding of the nature 
of interaction responsible for superconductivity (SC) in 
MgB-2- It has become clear that almost all facets of 
the phonon-mediated electron-electron interaction have 
a dramatic influence over the superconducting behavior 
of MgB2- For example, the electron-phonon matrix ele- 
ments change considerably as one moves away from the 
cylindrical Fermi sheets along r to A fj 16 1, anharmonic 
effects |6],|l5| have to be included in the dynamical matrix 
Jl6| , and finally k-dependent fully anisotropic Eliashbcrg 
equations have to be solved |L6| for a complete and ac- 
curate description of the superconducting properties of 
MgB 2 - Such a strong dependence of the superconducting 
properties of MgB2 on different aspects of the interaction 
has opened up the possibility of dramatically modifying 
its superconducting behavior by changing the interaction 
in various ways and thereby learning more about the in- 
teraction itself. Alloying MgB 2 with various elements 
and then studying their SC properties offers such an op- 
portunity. 

There have been several studies of changes in the SC 
properties of MgB 2 upon substitutions of various el- 
ements such as Be, Li, C, Al, Na, Zn, Zr,Fe, Co, Ni, 
and others ||,[l7] |2(]] . The main effects of alloying are 
seen to be (i) a decrease in transition temperature, T c , 
with increasing concentration of the alloying elements al- 
though the rate at which the T c changes depends on the 
element being substituted, (ii) a slight increase in the T c 
in case of Zn jL9| , pO[ substitution while for Si and Li 
the T c remains essentially the same, (iii) persistence of 
superconductivity up to x ~ 0.7 in Mgi- x Al x B 2 HfnJ, 
(iv) a change in crystal structure and the formation of a 
superstructure at x = 0.5 in Mg\_ x M x B 2 p7j , and (v) 
a change in the lattice parameters a and c. 

In an effort to understand the changes in the elec- 
tronic structure and the superconducting properties of 
MgB 2 alloys, we have carried out ab initio studies 
of Mgi- x Al x B 2 , Mgi- x Li x B 2 and Mg X - x Zn x B 2 al- 



loys. We have used Korringa-Kohn-Rostoker coherent- 
potential approximation |^l| , p2| in the atomic-sphere ap- 
proximation (KKR-ASA CPA) method for taking into ac- 
count the effects of disorder, Gaspari-Gyorffy formalism 
p3| for calculating the electron-phonon coupling constant 
A, and Allen-Dynes equation for calculating T c in 
Mgi^ x Al x B 2 , Mgi^ x Li x B 2 and Mg X - x Zn x B 2 alloys as 
a function of Al, Li and Zn concentrations, respectively. 
We have analyzed our results in terms of the changes in 
the spectral function p3] along r to A evaluated at the 
Fermi energy, Ep, and the total density of states (DOS), 
in particular the changes in the B p contribution to the 
total DOS, as a function of concentration x. 

For examining the possibility of superstructure forma- 
tion at x = 0.5, we have used ABINIT code Q, based 
on psuedopotentials and plane waves to optimize the 
cell parameters a and c as well as relax the cell-internal 
atomic positions of MgAlB^, MgLiB^ and MgZnB^ 
in P6/mmm structures. We have used these atomic 
positions to carry out a total energy comparison us- 
ing KKR-ASA CPA between the ordered and the sub- 
stitutionally disordered Mg\^ x Al x B2, M g\- x Li x B<i and 
Mg\- x Zn x B2 alloys at x = 0.5. Such an approach allows 
us to check the possibility of formation of a layered or a 
mixed superstructure at x — 0.5 in these alloys. Before 
we describe our results, we outline some of the computa- 
tional details. 

The charge self-consistent electronic structure of 
Mgi- x Al x B 2 , Mg x ^ x Li x B 2 and Mg\ .Zn.ll, alloys as 
a function of x has been calculated using the KKR-ASA 
CPA method. We have used the CPA rather than a rigid- 
band model because CPA has been found to reliably de- 
scribe the effects of disorder in metallic alloys pl| , p2| |. 
We parametrized the exchange-correlation potential as 
suggested by Perdew-Wang |^6| within the generalized 
gradient approximation. The Brillouin zone (BZ) inte- 
gration was carried out using 1215 fc— points in the ir- 
reducible part of the BZ. For DOS and spectral func- 
tion calculations, we added a small imaginary compo- 
nent of 1 mRy and 2 mRy, respectively, to the energy 
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FIG. 1. The calculated variation of T c as a function 
of concentration x in Mg\- x Al x B2, Mgi- x Li x B2 and 
Mgi- x Zn x B2 alloys. 

and used 4900 k-points in the irreducible part of the BZ. 
The lattice constants for Mg\_ x Al x B2, Mg\- x Li x B2 
and Mg\- x Zn x Bi alloys as a function of x were taken 
from experiments ]l7|-[l9[|. The Wigner- Seitz radii for 
Mg, Al and Zn were slightly larger than that of B. The 
sphere overlap which is crucial in ASA, was less than 10% 
and the maximum I used was l ma x — 3. 

The electron-phonon coupling constant A was cal- 
culated using Gaspari-Gyorffy Q formalism with 
the charge self-consistent potentials of Mg\- x Al x B2, 
Mg\- x Li x B2 and Mg\- x Zn x Bi obtained with the 
KKR-ASA CPA method. Subsequently, the variation of 
T c as a function of Al, Li and Zn concentrations was 
calculated using Allen-Dynes equation j24| . The average 
values of phonon frequencies u>i„ for MgB^ and AlBi 
were taken from Refs. iPJlcl respectively. For intermedi- 
ate concentrations, we took u>i n to be the concentration- 
weighted average of MgB2 and AIB2- For Mg\_ x Li x B2 
and Mgi- x Zn x B2 we used the same value of uji n as that 
for MgB 2 . 

The structural relaxation of MgALB^, MgLiB^ and 
MgZnB^ was carried out by the molecular dynam- 
ics program ABINIT with Broyden-Fletcher-Goldfarb- 
Shanno minimization technique [p5| using Troullier- 
Martins psuedopotentials [§?]], 512 Monkhorst-Pack Q 
k-points and Teter parameterization J25| for exchange- 
correlation. The kinetic energy cutoff for the plane waves 
was 110 Ry. 

Based on our calculations, described below, we find 
that in Mg\_ x Al x B2, Mg\_ x Li x B2 and Mg\- x Zn x B2 
alloys (i) the way T c changes depends on the location of 
the added/modified k-resolved states on the Fermi sur- 
face, (ii) the variation of T c as a function of concentration 



FIG. 2. The calculated spectral function along F to A, eval- 
uated at the Fermi energy, as a function of concentration x in 
Mgi- x Al x B2, Mgi- x Li x B2 and Mg\- x Zn x B2 alloys. Fig- 
ures (a) and (b) correspond to x = 0.1 and x = 0.3, respec- 
tively and the symbols open circle, open square, x and open 
triangle correspond to MgB2, Mg\- x Al x B2, Mg\- x Li x B2 
and Mg\- x Zn x B2 alloys respectively. In figure (c) the sym- 
bols open circle, open square and x correspond to AIB2, 
Mgo.1Alo.gB2 and Mgo.4Alo.eB2, respectively. For clarity, 
in figure (c) we have multiplied the spectral function of 
Mgo.4Alo.eB2 by 100. 

is dictated by the Bp contribution to the total DOS, (iii) 
at x — 0.5 a superstructure can form in Mg\_ x Al x B2 but 
not in Mg\^ x Li x B2 or Mg\^ x Zn x B2, and (iv) B layer 
shifts towards the impurity layer, more for Al than for 
Li or Zn. 

The main results of our calculations are shown in 
Fig. 1, where we have plotted the variation in T c of 
, M.,lh- Mgi- x Li x B 2 and Mg 1 - x Zn x B 2 alloys as 
a function of concentration x. The calculations were 
carried out as described earlier with the same value of 
fi* = 0.09 for all the concentrations. The T c for MgB 2 
is equal to ~ 30.8 K, which is consistent with the results 
of other works [0J^,|l^] with similar approximations. The 
corresponding A is equal to 0.69. For < x < 0.3, the 
T c increases slightly for Mg\^ x Li x B2 and M g\^. x Zn x B2, 
while it decreases substantially for Mg\_ x Al x B2, as is 
found experimentally [p|Jl7[]. Note that for x — 0.1, 
our calculation shows Mg\_ x Li x B2 to have a T c higher 
than that of Mg 1 - x Zn x B 2 by about 7K Jl|-|2^]. In 
Fig.l (right panel) we have shown the variation in T c 
in Mgi^ x Al x B2 alloys as a function of concentration for 
< x < 1. As a function of Al concentration, the T c 
decrease rapidly from 30 K at x = 0.05 to about 15 K at 
x = 0.1. The T c decreases slowly between x = 0.4 and 
x = 0.5. At x — 0.7 the T c vanishes and remains essen- 
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FIG. 3. The calculated total density of states at the Fermi 
energy (left panel) and the B p contribution to the to- 
tal DOS (right panel) in Mg\- x Al x B2, Mg\- x Li x B2 and 
Mgi- x Zn x B2 alloys. 

tially zero thereafter. The calculated variation in T c , as 
shown in Fig. 1, is in very good qualitative agreement 
with the experiments || ■ 

In order to understand the variation of T c in 
Mg 1 - X AI X B 2 , Mg 1 - x Li x B 2 and Mg 1 - x Zn x B 2 alloys as 
a function of concentration x, we have analyzed our re- 
sults in terms of the spectral functions, the contribu- 
tion of Boron p-electrons to the total DOS and the total 
DOS. In Fig. 2(a)-(c), we show the spectral functions 
along r to A direction evaluated at Ep in Mg\^ x Al x B 2 , 
Mg\_ x Li x B 2 and Mg\- x Zn x B 2 alloys for x = 0.1 (Fig. 
2(a)), x = 0.3 (Fig.2(b)), and x = 0.6 - 1.0 (Fig. 2(c)). 
From Figs. 2(a)-(b) it is clear that the substitution of Al 
in MgB 2 leads to creation of more new states along Y to 
A direction than the substitution of Zn or Li. Since the 
hole-like cylindrical Fermi sheet along V to A contributes 
much more to the electron-phonon coupling jl6| , the cre- 
ation of new electron states alongT to A direction weak- 
ens considerably the overall coupling constant A, which, 
in turn, reduces the T c more in Mg\- x Al x B 2 than in ei- 
ther Mg\- x Zn x B 2 or Mg\- x Li x B 2 . Thus, in our opin- 
ion, the way T c changes in MgB 2 upon alloying depends 
dramatically on the location of the added/modified k- 
resolved states on the Fermi surface. 

Having explained the differences in behavior of MgB 2 
upon alloying with Al, Li and Zn, we now try to un- 
derstand the changes in their properties as a function 
of concentration x. In Fig. 3(a)-(b) we have shown the 
total DOS at Ep (Fig. 3(a)) and the Bp contribution 
to the total DOS at E F (Fig. 3(b)) in M gi _ x Al x B 2 , 
Mg\- x Li x B 2 and M g\_ x Zn x B 2 alloys as a function of 
concentration x. We find that as a function of concentra- 
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FIG. 4. The calculated total density of states (solid line) 
for M go. 9AZ0.1-B2 (upper panel) , Mgo.gLio.1B2 (middle panel) 
and Mgo.9Zno.1B2 (lower panel) alloys. For comparison the 
total DOS for MgB2 (dashed line) is also shown. The dashed 
vertical line indicates the Fermi energy. 

tion, the variation in T c , as shown in Fig.l, follows closely 
the behavior of the total DOS at Ep and in particular 
the variation in Bp contribution to the total DOS at Ep. 
It is also not surprising to see that the vanishing of super- 
conductivity in Mgi- x Al x B 2 at x ~ 0.7 coincides with a 
very small Bp contribution to the total DOS. 

In Fig. 4(a)-(c) we show the total DOS of 
Mg x ^ x Al x B 2 , Mg x _ x Li x B 2 and Mg 1 - x Zn x B 2 alloys, 
respectively, at x = 0.1. In the same plot we also show 
the total DOS of MgB 2 obtained using the same ap- 
proach. The overall downward (upward) movement of 
the total DOS in Mgo.QAlo.1B2 (Mgo.gLi .iB 2 ) with re- 
spect to that of MgB 2 is due to the addition (removal) 
of electrons. In Fig. 4(c), the peak in the total DOS at 
around 0.53 Ry below Ep is due to the 3d states of Zn. 

Finally, we discuss the possibility of superstructure for- 
mation at x = 0.5. In Table I we show the optimized 
cell parameters as well as the cell-internal relaxations for 
MgAIBi, MgLiB/x and MgZnB± calculated using the 
ABINIT program. We find that the Boron layer shifts 
significantly more towards Al layer in MgAlB^ than to- 
wards either Li or Zn layer in MgLiB^ or MgZnB±, 
respectively. The shift of B layer by ~ 0.24 a. u. to- 
wards Al layer in MgAlB^ compares well with the corre- 
sponding shift obtained in Ref. (29). However, the shift 
of B layer towards Li and Zn layers in MgLiB^ and 
MgZnBi respectively implies that it is not simply due 
to the extra positive charge on the impurity layer, as 
suggested in Ref. |2{| in the case of MgAlB^. In Table I 
we have also listed the calculated ordering energy, E orc i, 
which is the difference between the total energies of the 



TABLE I. The calculated lattice constants a and c/a, the 
shift 8 of the B layer along the c-axis towards the impurity 
layer, and the ordering energy E or d at x = 0.5. The lattice 
constant a and the shift 5 are in atomic units while E or d is in 
raRy I atom. 
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c/a 
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MgAIBi 


5.799 


2.242 


0.24 


-12.1 


MgLiB 4 


5.685 


2.287 


0.04 


+4.8 


MgZnBi 


5.789 


2.254 


0.08 


+1.2 



ordered MgAlB^, MgLiB^ and MgZnB\ and the cor- 
responding disordered Mgo.5Alo.5B2, Mgo.5Lio.5B2 and 
Mgo.5Zno.5B2 alloys, obtained using KKR-ASA CPA 
method. It clearly shows the possibility of formation of a 
superstructure in Mgo.5Alo.5B2 because the fully-relaxed 
MgAlB^ is lower in energy by 12 mRy / atom in compari- 
son to the disordered Mgo.5Alo.5B2- However, within the 
limitations of our approach, we find that Mgo.5Lio.5B2 
and Mgo.5Zno.5B2 are unlikely to form superstructures 
since E or< i is positive in these two cases. Our results also 
show that a structure made up of layers consisting of a 
random mixing of Mg and Al atoms and described by 
CPA, is higher in energy than a structure made up of 
alternate layers of Mg and Al atoms . 

In conclusion, we have shown that in Mg\^ x Al x B2, 
Mg\- x Li x B2 and Mg\- x Zn x B2 alloys (i) the way T c 
changes depends on the location of the added/modified 
k-resolved states on the Fermi surface, (ii) the variation 
of T c as a function of concentration is dictated by the 
Bp contribution to the total DOS at Ep, (iii) at x = 0.5 
a superstructure can form in M g\- x Al x B2 but not in 
Mgi^ x Li x B2 or Mg\_ x Zn x B2, and (iv) B layer shifts 
towards the impurity layer, more for Al than for Li or 
Zn. 
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